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ABSTRACT 

Context. The external destruction of protoplanetary discs in a clustered environment acts mainly due to two mechanisms: gravitational 
drag by stellar encounters and evaporation by strong stellar winds and radiation. It is a fundamental question whether either of these 
mechanisms is important for the stellar evolution and planet formation process in young star clusters. 

Aims. We focus on the effect of stellar encounters in young dense clusters and investigate whether there are any observables that 
could trace this mechanism and its impact on disc evolution. If encounters play a role in disc destruction, one would expect that stars 
devoid of disc material would show unexpectedly high velocities as an outcome of close interactions. We want to quantify this effect 
by numerical simulations and compare it to observations. 

Methods. As a model cluster we chose the Orion Nebula Cluster (ONC). We reanalyzed observational data of the ONC to find en- 
counter signatures in the velocity distribution and a possible correlation with signatures of circumstellar discs. We use the nbody6++ 
code to model the dynamics of an ONC-like cluster and analyze the velocity distribution and the disc-mass loss due to encounters. 
Results. We found from the observational data that 8 to 18 stars leave the ONC with velocities several times the velocity dispersion. 
The majority of these high- velocity stars are young low-mass stars (t < 10 5 yr, m ~ 0.2-0.3 M Q ), among them several lacking infrared 
excess emission. Interestingly, the high-velocity stars are found only in two separate regions of the ONC - i) close to the cluster centre 
and ii) in the outer cluster region. Our simulations give an explanation for the location of the high-velocity stars and provide evidence 
for a strong correlation between location and disc destruction. 

Conclusions. The high-velocity stars can be explained as the outcome of close three-body encounters; the partial lack of disc sig- 
natures is attributed to gravitational interaction. The spatial distribution of the high-velocity stars reflects the initial structure and 
dynamics of the ONC. Our approach can be generalized to study the evolution of other young dense star clusters, like the Arches 
cluster, back in time. 

Key words, stellar dynamics - methods: N-body simulations, observational - stars: pre-main sequence, circumstellar matter 



1. Introduction 

According to current knowledge, planetary systems form from 
the accretion discs around young stars. These young stars are in 
most cases not isolated but are part of a cluster (e.g. lLada & Ladal 
2003). Densities in these cluster environments vary consider- 
ably, spanning a range of 10 pc~ 3 (e.g. 77 Chameleontis) to 
10 6 pc~ 3 (e.g. Arches Cluster) . It is still an open question as 
to how far interactions with the surrounding stars influence 
planet formation in dens e young clusters (« > 10 4 pc~ 3 ). These 
discs disperse over time (Haisch et al. l200ltlHii lenbrand 2002; 



Sicili a-Aguilar et all2006l ; ICurrie et alj|2008l) and in dense clus- 
ters the disc freq uency seems to be lower in the core (e.g. 
iBalog et all 120071) . This is attributed to external violent pro- 
cesses like photo-evaporation or encounter-induced disc mass 
loss, mainly caused by the massive stars that are concentrated 
in the cluster core. 

Earlier numerical investigations seemed to indicate that pho- 
toevapora tion should by far dominate the extern al disc de- 
struction dScallv&Clarkd l2001; Ada ms et alj|2004l) . The ques- 
tion of whether encounters play a vital role in the forma- 
tion process of stars and planets has be en studied far less 
and is still open (e.g. Ad ams et alj |2006). Only recently has 
it been shown that stellar encounters do have an effect on the 



discs s urrounding stars in a young dense cl uster (Olczak et al. 
1 2006 1 ; IPfalzner et al] 120061; IPfalznerl l2006t iMoeckel & Ballvl 
2006, l2007allbt IPfalzner & Olczakl l2007allbl; see also the re- 
view by IZinnecker & Yorkg 120071) . The importance of en- 
counters has been underestimated because in previous stud- 
ies th e focus has been on encounters betw een solar mass stars 
(e.g. IClarke&Prin"gl3ll991t lHellerlll995l) . However, discs are 
most affected when the masses of the sta r s involved in the 
encoun ter are unequal dOlczak et alj 120061 ; lMoeckel"&^ B ally 
2007b). Moreover, the massive stars in the center of a stellar 
cluster act as gravitational foci for the lower mass stars and are 
thus subject to repeated encounters. These accumulated pertur- 
bations can lead to a total destruction of a massive star's disc 
( IPfalzner et al.ll2006l) . 

The main focus of this paper is the question of whether there 
is direct observational evidence for encounters among star-disc 
systems in young clusters and for disc-mass loss due to encoun- 
ters. It is difficult to distinguish observationally whether photo- 
evaporation or gravitational interaction are responsible for the 
loss of (outer) discs. The reason is that in both cases interaction 
with a massive star is the most destructive process. Thus the ob- 
servation of a decreased disc frequency in the cluster core, as 
mentioned above, does not allow one to favor either of the two 
mechanisms. The difficulty of tracing stellar encounters directly 
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by observations is their short duration. Direct imaging of tidal 
tails of a disc would be the best proof but the probability of such 
an event is very low: the perturbed disc circularizes and tidal tails 
dissipate on a time scale of <1000yr. Nevertheless, there ex- 
ist ob servations of tida l tails or spiral arms in star-disc sy stems 
(e.g. iBeust etal.ll200l ICabrit et alj|2006t iLin et alj 12000) . but 
often it is unclear whether these can be attributed to the passage 
of an unbound perturber or a binary companion or gravitational 
instabilities caused by a giant planet. Due to the short dissipa- 
tion time scale of encounter-induced tidal tails, there would be a 
high probability for the encounter partner still to be located close 
(< 10 4 AU) to the perturbed disc. However, source confusion in 
dense young clusters like the ONC, where stellar encounters are 
frequent enough to study the effect on stellar discs, would pre- 
vent a clear identification. Set in relation to the rate of encounters 
in which prominent tidal tails may be formed, we estimate an up- 
per limit of just four candidate stars for a direct observation of 
tidal tails in a protoplanetary disc due to encounters in the ONC 
(see AppendixlAlfor a detailed calculation). However, an unam- 
biguous imprint of an encounter among stars is the high velocity 
of a star which has been expelled in a close gravitation al inter- 
action , mostly as a result of a three-body encounter (see Heggie 
fl975l) . Thus an analysis of the velocity distribution of a cluster is 
the key to finding candidates of close encounters between young 
stars. 

Here the ONC is used as a prototype young cluster because 
it is one of the best-studied regions in our galaxy, and the only 
young dense cluster for which velocit ies of more than 10 00 of 
its members have been determined dJones & Walked Il988h . In 
addition, its high density suggests that stellar encounters might 
be relevant for the evolution of circumstellar discs. Throughout 
this work we assume that initially all stars are surrounded by 
protoplanetary discs. This is justified by observations that re- 
veal dis c fractions of nearly 100% in very young star clus- 
ters (e.g.lHaisch et al.l2000l:lLada et al.l2000tlHaisch et all200lt 
lHillenbrandll2005l) . 

In §|2] we present results from a search for candidate stars 
of close encounters in the publicly available observational data 
of the ONC . For this purpo s e we h ave reanalyzed the inves- 
tigations of Uones & Walker d!988l) . iHillenbrandl d!997l) . and 
iHillenbrand et al.l 0998) for stellar velocities and infrared ex- 
cess. The basic properties of the ONC used for our numerical 
model are described in |3] There we also discuss theoretical es- 
timates and observational data related to binary populations and 
encounters. Afterwards we present results from a numerical ap- 
proach to this problem (@. Observational and numerical results 
are summerized and discussed in $5] 

2. Observational data 

Our search for candidate stars of close encounters from obser- 
vational data concentrates on velocity surveys of Orion. There 
exists only one large data set of the ONC that provides velocity 
in formation for several h undred stars, the proper motion survey 
ofQ ones & Walker! d 19881). Fortunatel y, the same stars have been 
analyzed by IHillenbrand et all d!998l) for infrared excess emis- 
sion as an indicator of circumstellar material. They defined a 
quantity measuring the magnitude of the near-infrared excess, 
A(/ c - K), 



A(/ C - K) = (I C - ^observed - 0.5A V - (I C ~ K) 



photosphere •> 



(1) 



where the first term is the observed Iq-K color, the second term 
the contribution of reddening calculated from the extinction val- 
ues derived from V — Iq colors as in Hillenbrand (1997), and the 



third term the contribution of the underlying stellar photosphere. 
This measure of the near-infrared excess, A(/c - K), traces only 
the innermost (<0.1 AU) part of the circumstellar disc. Hence, 
the absence of such emission is not to be confused with a com- 
plete absence of a protoplanetary disc, yet is dependent on a 
number of parameters, such as disc a ccretion rate, inclination , 
inner hole, and stellar mass and radius (Hillenbr and et al.f i 998). 
However, in the following, we will refer to stars that lack infrared 
excess emission as "discless", emphasizing that these stars may 
be well surrounded by circumstellar material but do not show the 
typical infrared excess of a young star-disc system. 

In addition Hillenbrand ( 1997) has investigated stellar prop- 
erties like mass and age. The investigation is based on opti- 
cal photometric and spectroscopic data and covers only about 
half of the stellar population of the ONC, while the more em- 
bedded stars are n ot accessible at this wavelength. However, 
IHillenbrand! ([1997) states that the investigated stellar sample is 
representative of the entire population of the ONC. Masses and 
ages were derived from an HR diagram via comparison with 
theoretical pre-main sequence evolutionary tracks. This method 
leaves some uncertainty as to the absolute stellar age (~0.5 dex, 
see Appendix IB1 and mass calibrations, with mass deviations of 
about factors of two and age differences of several Myr among 
different models. Additionally, photometric errors translate into 
uncertainties of the derived masses, but this is significant only 
for stars more massive than M«1.5M . However, all these 
sources of uncertainties are less relevant for the present inves- 
tigation because (i) the stars of interest are of low mass as will 
be shown later, and (ii) the derived conclusions are based mainly 
on relative masses and ages. 

We have merged the data of the three investigations de- 
scribed and excluded stars for which measurements of proper 
motion or infrared excess are missing, resulting in a database of 
655 stars. In order to achieve the most secure distinction possible 
between stars with infrared exc ess and pure photospheric emis- 
sion, we adopted the criterion o flSicilia-Aguilar et alj d2005l) and 
excluded all stars that have an infrared excess A(/c - K) in the 
range < A(/c - K) < 0.5. A star is classified as discless if 
A(/ c - K) < 0, or as a star-disc system if A(/ c - K) > 0.5. 
This additional selection criterion reduces the sample of stars 
that have been used for this investigation to 405, among them 
266 star-disc systems and 139 stars lacking any excess emission. 

In Fig. [1] we show the velocity distribution of stars with and 
without discs, adopting a distan ce of 420 pc to the ONC (cf. 
lJeffriesll2007l;lMenten et alj|2007l) . The velocity distribution has 
been built by binning each component of the proper motion, /u x 
and fly, of each star separately, not its total two-dimensional mo- 
tion, or in other words by summing the Gaussian velocity distri- 
bution of each spacial direction, x and y, which again results in 
a Gaussian. The reason not to bin the total two-dimensional mo- 
tion is the improvement of statistics due to the twice as large data 
set. For comparison, a Gaussian with a one-dime nsional velocity 
disper sion <x = crj* = 2.5 km s _I , as derived by iJones & Walker! 
(1988), is superimposed (dashed line). The distribution shows 
the expected features: At velocities fi xy < 3o-™ the shape is ap- 
proximately Gaussian, though a distinct peak at zero velocity is 
present. This peak is only prominent for the x-components of the 
proper motion but is independent of the applied binning. We sup- 
po se that this is an a r tifact of the plate reduction technique used 
by IJones & Walker! d!988h . At higher velocities, fx x , y > 3cr J *, 
there is an overabundance of stars when compared with the theo- 
retically expected Gaussian distribution. In this regime stars are 
not bound energetically to the cluster (see Bin nev & Tremaina 
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Fig. 1 . Velocity distribution from observational data, adopting a distance of 420 pc to the ONC (cf. I Jeffries! 120071: iMenten et al.l 
l2007h . See text for details of its construction. The bin width is 0.7 km s . The stars have been divided into two groups according 
to the disc signature A(/c - K) (see text): grey boxes represent star-d isc systems, lines represent discless stars. For comparison a 



Gaussian with dispersion <x = crj* = 2.5 km s 1 ( Jones & Walker 1988) is superposed (dashed line), a) Complete samp le of stars in 



the ONC with available proper motions from lJones & Walker! d 198 8) and infrared excess from lHillenbrand et alj (QJ98). b) Like a) 
but zoomed into the lower part, c) Like b), but with restricted stellar ages (see text), d) Like b), but with restricted stellar ages and 
velocity errors as described in the text. 



119871 Eq. 8-3). In the following we will refer to these stars as 
"high-velocity stars". 

The number of high-velocity stars, M, V s,obs = 50, is remark- 
ably high, much higher than the theoretically expected A^vs 5 15, 
as derived in § 13.21 In fact, A^vs, obs should be even lower than 
the theoretical A^vs because only about 1 /5 of the ONC stellar 
population is covered by the observational data. The reason is 
most likely a contribution from foreground stars, which show 
large proper motions due to projection effects. In order to ex- 
clude probable foreground stars from the sample of ONC stars, 
we plot in Fig.|2]age vs. velocity for the high-velocity stars, again 
for the two groups of discless stars and star-disc systems. Two 
opposite trends are visible: The discless stars with ages t > 1 Myr 
have a ~10 times higher maximum velocity than younger ones 
(~200kms _1 compared to ~20kms _1 ) and a twice as high mean 
velocity (~13kms _1 compared to ~30kms~'). This strong cor- 
relation is likely due to a large fraction of foreground objects 
among the stars with ages t > 1 Myr, and since this age boundary 
is consistent with the mean cluster age, it marks a conservative 
upper limit to the age of probable cluster members. Conversely, 
the star-disc systems are similarly distributed in both age groups, 
with comparable maximum and median velocities (~25kms _1 
and ~llkms _1 , respectively). However, we exclude all stars 
with ages t > 1 Myr in the further discussion for three reasons: 



(i) the age of these stars exceeds our simulated cluster age, so 
they are dynamically distinct, (ii) at least some of these stars 
might not be members of the ONC due to their high age (and 
velocity), and (iii) protoplanetary discs with ages t > 1 Myr may 
be subject to significant evolution due to internal processes 1 , so 
that the effect of encounters could not be quantified. 

Of the group of high-velocity stars younger than 1 Myr, the 
"probable high-velocity cluster members", about half show disc 
signatures, the rest does not. The fact that stars with velocities 
of several tens of kms -1 are not present in this group is in ac- 
cordance with dynamical estimates: the probability of both the 
generation and detection of stars with such high velocities is very 
low due to the need for a very close approach and short traverse 
of the cluster. 

The uncertainties in the estimated ages of stars with and 
without infrared excess emission do not alter the fact that some 
of the youngest high-velocity stars lack infrared excess. This is 



1 According to studies of p roperties and evolu t ion of protoplanetary 
discs in young clusters (e.g. Hillenbrand 2005; Sicilia-Aguilar et al. 
2006), it is valid to assume that the protoplanetary discs in the ONC 
have been only marginally subject to internal disc processes. As such, 
we assume in the following that the effect of external processes on discs, 
i.e. photevaporation and encounters, have not been masked by internal 
processes at the current age of the ONC. 
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Fig. 2. Distribution of age and velocity of the high-velocity stars 
in the ONC, classified as star-disc systems (open circles) and dis- 
cless stars (filled circles). The mean age of the ONC, t — 1 Myr, 
and the minimum one-dimensional velocity of the high-velocity 
sample, 3<x^, are indicated by the horizontal and vertical dashed 
line, respectively. The uncertainty of the determined ages of 
~0.5 dex (see Appendix IE1 for the derivation) is indicated by the 
vertical error bar. At high velocities the (horizontal) velocity er- 
ror bars are smaller than the symbol size. 



noteworthy because in the early stages of stellar evolution one 
would expect accretion rates to be high and circumstellar mate- 
rial to be close enough to the star that significant excess emission 
can be detected. However, because the excess emission mea- 
sured by A(/c - K) originates close to the stellar surface, it is 
strongly dependent on the geometry and orientation of the disc. 
A more robust a nd sensitive indicator of circumstella r discs is 
the/: 



L color dMever et al.lll997tlHaisch et al.l 2000), tracing 

material out to radii of ~0. 1 AU (IHaisch et al. f l2005l) . We have 
thus cross-checked the "discless" stars with A(/c - K) < for 
excess emission at longer wavelengths. The results are shown in 
Table [U 

Two of the stars previously determined as "discless" do show 
a typical emission signature of warm circumstellar matter. Four 
high-velocity stars remain that lack infrared excess. We cannot 
determine whether the pure photospheric colors point to the ab- 
sence of a circumstellar disc. We will discuss this possibility 
later. 

We have plotted the positions and velocity vectors of the se- 
lected probable high-velocity cluster members in Fig. [3^. Two 
features are apparent: (i) Most high-velocity stars are concen- 
trated in the inner tenths of a parsec around the most mas- 
sive ONC member, 8 l C Ori, and (ii) "outliers" are located more 
than 1 pc from the cluster center with some stars moving in 
radial directions away from the cluster centre. Although one 
would expect the former, as most encounters usually happen in 
the dense cluster centre, there are two unexpected features: (i) 
The region between the cluster centre and outer cluster parts, 
0.4pc<r< 1.0 pc, is devoid of high-velocity stars, and (ii) the 
velocity vectors of several stars in the outer cluster parts do 
not point away from the cluster centre (as if these stars had an 
encounter in the less dense cluster parts). In order to be confi- 
dent that the encounter candidate stars are indeed high-velocity 
stars and not affected by observational uncertainties, we restrict 
our sample of high-velocity stars to maximum velocity errors 
of less than 30% or below 2.0kms~'. This reduces the num- 
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Fig. 3. Positions and velocities of high-velocity star-disc systems 
(open circles) and discless stars (filled circles) in the ONC. The 
frames are centered on the most massive star, l C Ori, which is 
marked by a cross. Thin lines indicate the distance (in pc) a star 
would have moved in a period of 2 • 10 4 yr. Thick lines and the 
opening angle reflect the error of the magnitude and the direc- 
tion of the proper motion, respectively, a): Sample of all high- 
velocity stars with ages < 1 Myr. b): Restricted sample of high- 
velocity stars with ages < 1 Myr: only stars fulfilling the con- 
straints on the velocity errors as dicussed in the text are shown. 



ber of high-velocity stars to a total of only eight, among them 
just one star without infrared excess, JW 45. The correspond- 
ing position-velocity diagram is shown in Fig. [3}) (see also the 
corresponding velocity distribution in Fig. [TJl). However, apart 
from the concentration of high-velocity stars close to the cluster 
centre, the same features as in Fig. [3^ are apparent. The reason 
for the different central concentration is the strong acceleration 
of stars close to the cluster centre which results in large errors 
of the derived proper motions. Consequently our velocity error 
criterion preferentially excludes stars close to the cluster centre. 
We will address the particular features in §0 where we present 
the results from numerical simulations of our dynamical model 
of the ONC. 
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" Stellar ID from lJones & Walked i!988h 

Projected distance from 6 l C Ori. 
c Proper motion and associated error. 
d Probability of cluster membership. 
' Apparent magnitude in specified band filter. 
' Spectral type. 
s Stellar mass. 
* Stellar age. 

' Detection of excess emission. 
Table 1. Properties of discless high-velocity stars in the ONC. 



Below we review the known properties of the ONC, which 
we use to model the cluster in an as realistic a way as possible. 

3. Structure and dynamics of the ONC 

The ONC is a rich stellar cluster with about 4000 mem- 
bers with masses ot>0 .08 Mq in a volume ~5pc across 
dHillenbrand & Hartmannl 119981: IHiUenbrand & Carpenter] 
2000). Most of the objects are T Tauri stars. The mean stellar 
mass is about in » 0.5 M Q and the half-mass radius Rh m w 1 pc 
(IScall y et al. 2005). Recent studies of t h e stellar mass distri- 
bution dmienbrand & Carpenterl [20001: iLuhman et all 120001: 
iMuench et all 120021: ISlesnick et alj|2004l) reveal no significant 



deviation from the generalized IMF of IRroupal (12002). 



-1.3 



0.08 < m/M < 0.50, 
O.5O<m/M < 1.00, 
1.00 < m/M < oo. 



(2) 



The shape of the system is not perfectly spherical but elon- 
gated in the north-south direction. The probable reason for 
this asymmetry is the gravitational potential of a massive 
molecular ridge in the backg round of the cluster, OMC 1 
(Hillenbr and~& Hartmannl fl998). The mean age of the whole 
cluster has been estimated to be frjNC ~ 1 Myr, alt hough a sig- 
nificant age spread of the in dividual stars is evident (Hillenbrand 
ll997t|Palla & Stahlerl2000l) . Radio observations bv lWilson et afl 
( 1 19971) show that only a few solar masses of ionized gas are 
present in the inner ~ 1 pc. 

The density and velocity distribution of the ONC resembles 
an isothermal sphere. The central n umber density p C ore in the 
inner 0.053 pc reaches 4.7 • 10 4 pc" 3 dMcCaughrean et al.ll2002h 
and makes the ONC the densest nearby (< 1 kpc) young stel- 
lar cluster. The dense inner part of the ONC, also known as 
the Trapezium Cluster (TC), is characterized by /?xc^0.3pc 
and A^rc ~ 750 , or »tc ~ 10 3 PC" 3 . In their proper motion study 
of the ONC, iJones & Walked dl988l) found the velocity dis- 
persion to be nearly constant at all cluster radii and ob- 
tained a one-dimensional velocity dispersion crj^ = 2.5kms~ 1 . 
This translates into a three-dimensional velocity dispersion of 
cr™= V3cr j w= 4 31^^1 Recently, iFuresz et all d2007l) ob- 
tained a somewhat higher one-dimensional velocity disper- 
sion of cr^p = 3.1 kms -1 from radial velocity measurements. 



However, they caution that their velocity distribution of ONC 
stars has a peak that is too low compared to the expected 
Gaussian di stribution with dispers ion cr^, so we will rely on 
the result of lJones & Walked" (119881) . 

Using this value, the virial ratio Q vu of the ONC becomes 



Gvir = 



1 ^°3D ) 

~2GM 



1.5, 



(3) 



where M = ;«A^~2OOOM . This indicates that the ONC seems 
to be gravitationally unbound (Q vlr > 1). However, the gas mass 
of th e background molecu lar cloud OMC 1 should be consid- 
ered. |Huff&JStahlej] (|200a) suggest that the molecular cloud that 
formed the ONC contained about 6000 M . The remaining part 
of it in the cluster backg round still has a mass of about 2000 M Q 
dHillenbrand et al.lll998h . 

In the most re cent study on ci rcumstellar discs in the 
Trapezium Cluster, lLada et a l. (2000) found a fraction of 80- 
85 % discs among the stellar population from the L-band excess. 
This is in agree ment with an earlier inve stigation of the complete 
ONC in which IHiUenbrand et all (Q998) report a disc fraction of 
50-90 % (though relying only on Iq - K colors) and justifies the 
assumption of a 100 % primordial disc fraction. 

The treatment of binaries in our model of the ONC is cru- 
cial, since they can have a strong impact on the evolution of 
cluster dynamics. In the next section we give a summary of the 
main properties of the ONC binary population that determine 
our model setup. 

3.1. Binaries in the ONC 

In the ONC the binary rate for solar-type stars is ~50%. From 
observations alone we have only a very limited knowledge of 
the distribution of binary periods, eccentricities or mass ratios of 
the ONC. However, combining observational data and numer- 
ical simulations, the initial properties of the primordial binary 
populat ion in a stellar aggregate have be en modelled bv lKroupal 
( 1995b} and Eouwenhoven etal] d2007l) . These can be applied 
to some degree to the ONC. 

The investigation of IRroupal (Il995bb is based on the prop- 
erties of the Taurus-Auriga binary population and constructs the 
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primor dial population b y inverse dynamical population synthe- 
sis (see iKroupa 1995a) and pre-main-sequence eigenevolution. 
The resulting distribution s are approximately th e log-n ormal pe- 
riod distribution fp(P) of Duq uennov & MavoJ d!991l) . 



f P (P) oc exp 



log P - p P 



2a% 



<P<P n 



with mean pp = log P- 4.8 and standard deviation 
(7> = c\ og p — 2.3, P in days, a thermally relaxed eccentric- 
ity distribution f e (e), 



f e (e) = 2e, 



< e < 1, 



(5) 



and a mass ratio distribution f q (q) obtained by random pairing 
of stars, 



fqia) k q 



7q 



(6) 



where q = Mi/Mx, M\ the primary, M2 the secondary mass, and 
y q - a, a the slope of the mass function of the stellar system. 

The log-normal period distribution fp(P) results in an ap- 
proximately log-normal semi-major axis distribution /«(«), the 
shape of whic h is slightly dependent on the distribution over bi- 
nary mass M jKouwenhoven et aiT 2007): 



log a = -logP 



3 B \2GM) 



°"ioga - ^o- iog p. (7) 



iKouwenhoven et all d2007l) have analyzed the current binary 
population of Scorpius OB2, under the reasonable assumption 
that it is still close to its primordial state. Accounting for differ- 
ent observational biases by means of comparison with simulated 
observations of model associations, they recovered a somewhat 
different primordial binary population: The semi-major axis dis- 
tribution of Sco OB2 is most consistent with a fiat distribution in 
logarithmic space, and is equivalent to 



f a (a) oc a 7 



< a < a„ 



(8) 



with Ami,, s; 5 R , flma x ~5 • 10 6 R q, and y a = —\, which is also 
known as Opiks law dOpikll 1924ft . The eccentricity distribution 
could not be well constrained, but the observations are consis- 
tent wit h a thermal distribution, given by Eq. ((5]). Unlike Kroupa 
dl995bl) . IKouwenhoven et ail d2007l) find a power law depen- 
dence of the mass ratio distribution with y q m -0.4. This is much 
flatter than Eq. (O and favors massive companions for massive 
stars. 

The model of IKouwenhoven et al.l d2007l) seems more appli- 
cable to our case, as observational studies of the ONC favor a 
flat distribution of th e semi-major axes (e.g. lPadgett et all 19971; 
Reip urth et al.ll2.Q07h . Theoretical cons iderations bas ed on three- 
body encounters give similar results (Valto nenll 1997b . 

A thermal eccentricity distribution is expected from energy 
equip artition as a result of multiple soft encounters dHeggid 
1975) and is also found from observations, though only for bina- 
ries with separations a > 10-50 AU. Very close systems are sub- 
ject to circ ularization due to tidal effects occurring dur ing stellar 
evolution dDuquennov & Mavoilll 99 U iMathieull 19941) . 

The shape of the observed mass ratio distribution is not 
well constrained by observations. However, to a good ap- 
proximation the mass ratio distribution can be described by 
a p ower law as given by Eq. ^ over a wide mass range 
(e.g.llrimb le 1990; Mal kov & Zirmeckeiil200lt lRucinskill200U 
IShatskv & Tokovininl [20021: IValtonenl 120041) . We thus favor a 
single mass ratio distribution for primordial binaries g i ven b y 
Eq. (O with y q = -0.4 as derived bv lKouwenhoven et al.l d2007l) . 



Binaries have strong effects on the overall cluster dynamics 
mainly through close interactions with single stars or other mul- 
tiples. In the following we estimate the typical encounters that 
could generate high- velocity stars. 



(4) 3 2 Three-body encounters in the ONC 



As is well known, the non-hierarchical motion of three bodies, 
known as the three-body problem, has no analytical solution, and 
the chaotic motion of the members can only be investigated nu - 
merically in a statistical manner (see Valtonen & Mikkola 1991). 
In her study of triple systems with negative total energy. lAnosoval 
dl986) found that about 95 % of three-body systems decay after a 
close triple approach of the components. In most cases (~80 %), 
ejection leads to escape, but can also result in the formation of 
a hierarchical triple system, with one body in an extended orbit. 
The lowest mass member has the highest probability of being 
ejected, about 80%. 

In the following we assume that high-velocity stars are typ- 
ically generated in triple systems with negative total energy. We 
will justify this assumption later in this section. If the motion 
has not been significantly perturbed since the encounter, one 
can draw conclusions about the underlying encounter parame- 
ters from the dynamics of the ejected body. In the case of the 
ONC, it is valid to assume that high-velocity stars with veloci- 
ties v > 3cr^* « 13 km s _1 are effectively unperturbed before they 
escape from the cluster (see AppendixO. Then the time to reach 
the outskirts of the cluster is t esc ^R/v< 0.2 Myr. 

For a crude estimate of the compactness of the three-body 
system (with negative total energy) from which a member is 
ejected with f>3cr ^f, we evaluate the sca ling of the median 
escape speed from Sterzik & Durise3 d 19951) . 



<Vesc) ~ ~ 



1 / |£bl 



2 \<m eS c)/' 



(9) 



where (m esc ) denotes a weighted mean of the escaped particle 
masses, \Eo\ocMf ot /R is the total system energy, and M tot and R 
are the total mass and the scale length of the system. 

We assume that the encounter occurred in the dense 
Trapezium Cluster where it is most probable (see Fig. 4 of 
Pfalzner et al. 2006j). Due to mass segregation of the cluster, 
the mass of the most massive component of the three-body sys- 
tem is likely to be several times the mean stellar mass in the 
ONC; we adopt M to t =4M Q for the system mass. The mass of 
the ejected body is assumed to be half the mean stellar mass, 
(wiesc) = 0.25 M Q . With these assumptions we obtain R < 100 AU 
for the scale length of the system. However, the minimum ap- 
proach which eventually leads to an ejection will be much closer. 
Since the disc radius of a low-mass star is about rj * 100 AU, the 
ejected co mponent can lose m ore than 90 % of the disc mass (see 
Table 3 of lOlczak et aD2006l) . 

Initially we assumed that triple systems that generate high- 
velocity stars have negative total energy. In AppendixlDlwe show 
that this is a valid assumption if the ejected high-velocity stars 
do not exceed velocities of a few tens of kms . In fact, all of 
the high-velocity stars from numerical simulations and observa- 
tional data do not show higher velocities. A detailed analysis of 
our numerical simulations shows that in all cases the encounter- 
generated high-velocity stars are the lowest mass component 
ejected from a three-body system with a massive tight binary 
(Mbm ^ 20 M , flbin SS 50 AU), leaving on a nearly parabolic or- 
bit relative to one of the binary components. We thus conclude 
that the triple systems which generated the observed low-mass 
high-velocity stars must have been bound as well. 
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How many such encounters do we expect in the ONC? 
For simplicity we assume the three-body interaction to result 
from a single-binary encounter, where the binary has a semi- 
major axis a<100AU. As a volume relevant for close en- 
counters we consider the Trapezium Cluster. Referring to Eq. 
( IC.U the time scale for a three-body encounter at 100AU is 
?enc,iooAU ~ lOMyr. Since the time the star remains in the clus- 
ter volume is f esc 5 0.2 Myr, the probability of the detection of 
a high-velocity star is P bvs « * esc /*enc,iooAU £ 0.02. Knowing the 
number of stars located in the Trapezium Cluster, Njc, we ex- 
pect A^hvs ~ PhvsNjc ^ 15 high-velocity stars in the ONC. 

In summary, high-velocity stars with velocities v>3cr™ (or 
velocity components v XyV , > 3crjp ) have most probably experi- 
enced exactly one close (< 100 AU), disruptive encounter in the 
cluster centre and leave the cluster on a radial trajectory, i.e. with 
the velocity vector pointing away from the cluster centre, pre- 
serving the dynamical information of the encounter. We expect 
about one dozen such stars in the ONC. 

In the following we will compare the observationoal data and 
theoretical estimates with numerical simulations of a dynamical 
model of the ONC. 



4. Numerical simulations 

Th e basic dynamical m odel of the ONC used here is described 



Olczak et al.l (12006). with several extensions discussed 



iPfalzner & Olczakl d2007bl ). We summarize the main properties 
of the initial stellar distribution: the simulations start in virial 
equilibrium, Q = 0.5, with a radial density profile, p ct r 2 , a cen- 
tral density ~4-10 5 pc~ 3 , and a Maxwellian velocity distribution. 
Here we include additional effects like a varying background 
potential and a population of primordial binaries. We have per- 
formed 20 runs with different random configurations of positions 
and velocities from the given distributions to establish a statis- 
tically robust database. Unless explicitly declared otherwise the 
presented results refer to the whole set of runs. The reason to 
extend our basic model of the ONC is that we want to compare 
the velocity distribution to the observational data. Both the back- 
ground potential and primordial binaries have a large impact on 
the resulting velocity distribution. However, for comparison and 
analysis of the contribution of the background potential and pri- 
mordial binaries, we also performed simulations of a single par- 
ticle model, a single particle model with a background potential, 
and a cluster with a primordial binary population but without a 
background potential. 

A background potential increases the cluster virial mass and 
thus results in a higher velocity dispersion (cf. Eq. (O). In 
our simulations the background potential is represented by a 
Plummer sphere, 



P(r) 



3M 



1 



4na 3 (l + 4)5/2' 



(10) 



with mass M and characteristic length scale a. The Plummer 
sphere is set up with 6000 M initially and a mass loss rate 
M = 4OOOM Myr~ 1 . Due to the continuous mass loss, after 
1 Myr, 2000 M of gas is left. The exact time-scale and time 
dependency of the mass loss is not crucial, as long as the gas 
expulsion time scale is of the order of the dynamical time scale 
of the system, f exp « ti yn , which appears to be the case for the 
ONC (cf. §0. The length scale is set to a = 0.6pc to match the 
observed velocity dispersion at a simulation time of f s i m = 1 Myr. 

The effect of primordial binaries is more complicated. The 
interactions of binaries and single stars or other binaries have 
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Fig. 4. Three-dimensional velocity dispersion as a function of 
time for the previous (dashed) and new (solid) numerical mod- 
els of the ONC. The dotted horizont al line marks th e three - 
dimensional velocity dispersion from |Jones&JValker (1988), 
crjj* =4.3kms~', the vertical dotted line marks the mean age of 
the ONC, ?onc = 1 Myr. 



the potential to change the velocity distribution much more than 
a single star model. In particular, three body encounters between 
singles and binaries usually lead to the expulsion of the lowest 
mass member from the small A^-body system with a high veloc- 
ity. This mechanism is especially important for the high-velocity 
fraction of particles of interest in this study (see !33.2I >. 

We have se t up a series of simulations with the models of 
iKroupd d!995bl) and lKouwenhoven et ail d2007l) and found that 
results do not depend on the choice of one particular model. The 
initial binary frequency was chosen to be 75 %. Although the 
observed binary frequency in the ONC is ~50 %, it is necessary 
to start with a higher binary rate as about one third of the binaries 
are distroyed within the first 1 Myr due to dynamical evolution. 

In terms of global cluster dynamics, the new model with a 
background potential and primordial binaries provides a much 
better fit to the ONC data than the previous single star model 
without background potential. We illustrate this by means of 
the time-evolution of the three-dimensional velocity dispersion 
of both models in Fig. [4] In the previous model there was not 
enough mass confined in the cluster to reproduce the observed 



velocity dispersion of the ONC, cr^ w ■ 



4.3 kms , at ~1 Myr. 



The new model gives a much better result. The continuous, steep 
falloff after ~0.2 Myr is due to the response of the stellar system 
to the gas expulsion. In the following, we will discuss results of 
numerical simulations based on the new dynamical model of the 
ONC. 

Olczak et alj d2006l) . IPfalzner et all 



As demonstrated in 



(2006), and IPfalzner & Olczakl d2007al) . stellar encounters in 



dense clusters can lead to significant transport of mass and an- 
gular momentum in protoplanetary di scs. In the pre s ent inv es- 
tigation we have used Eq. (1) from Pfalzner et alj d2006l) to 
keep track of the disc-mass loss of each star due to encoun- 
ters. Our estimate of the accumulated disc-mass loss is an upper 
limit because the underlying formula is only valid for co-planar, 
prograde encounters, which are the most perturbing. A simpli- 
fied prescription assigns stars into one of two distinct groups: 
if the relative disc-mass loss exceeds 90% of the initial disc 
mass, stars are marked as "discless"; otherwise they are termed 
"star-disc systems". This approach is justified by the interplay 
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Fig. 5. Velocity distribution of cluster stars from simulations. 
The construction of the distribution is described in the text. The 
bin width is 0.7 km s -1 . The stars have been divided into two 
groups according to their extant disc mass (see text): boxes rep- 
resent star-disc systems, impulses represent discless stars. For 



comparison a Gaussian with dispersion <x = 
superposed (dashed line). 



ID 



;2.5kms is 



of three effects: (i) A disc-mass loss of this order lowers the 
density in the disc significantly, in particular in the outer parts; 
the disc size decreases, (ii) The accompanying angular momen- 
tum loss enhan ces accretion of the extant material onto the star 
(Pfalzner 2006). This leads temporarily to an increase of the in- 
frared excess but soon fades after a short intense accretion phase 
(Pfalzner (2008), subm. to A&A). (iii) The loose distribution of 
circumstellar matter lowers the shielding of the disc midplane 
against photoevaporation. The interplay of these effects leads to 
a fast dispersal of the disc material. From the observational point 
of view, the corresponding star would show pure photospheric 
emission on the order of some 10 3 yr after the encounter. 

In Fig. [5] the velocity distribution of cluster stars is shown 
after 1 Myr of evolution. This is done separately for discless 
stars and star-disc systems. For comparison with the observa- 
tional data presented in Fig. [lj, separate velocity distributions 
of two spacial directions (here x and y) have been added to 
mimic the distribution of proper motion data. Unless in wide 
sy stems, primaries and se condaries could not have been resolved 
bv lJones & Walkerl (1 19881) who worked with seeing-limited im- 
ages. So the presence of binaries in the numerical simulations re- 
quires a special treatment of velocities. Accounting for unequal 
mass components and nebulosity, we adopt a minimum separa- 
tion o f 1000 AU (corresponding to ~2.5"; see Mayne & Navlor 
2008, and references therein) for the visual resolution of a bi- 
nary system. For closer systems, only the primary component is 
taken into account and the center-of-mass velocity of the system 
is assigned. For wider systems, both components are treated as 
single stars. This prescription is simple and rough but appropri- 
ate to avoid the inclusion of large velocity components in tight 
binaries. 

The velocity distribution shows common features with 
Fig- 0J, which was obtained from observational data: the bulk 
of the stars forms a relaxed system which manifests in the ap- 
proximate Gaussian velocity distribution with a characteristic 
velocity dispersion cr = a J ^. Moreover, a small fraction of stars 
exists with much higher velocities v v . v > 3cr[p, the previously 
described "high-velocity stars". 
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Fig. 6. Positions and velocities of high-velocity star-disc systems 
(open circles) and discless stars (filled circles) from cluster sim- 
ulations. The frame is centered on the cluster centre. The lines 
indicate the distance (in pc) a star would have moved in a period 
of2-10 4 yr. 



The velocity distributions of stars that have lost their disc 
due to close encounters and those that have retained their disc 
differ. The discless population consists of a larger fraction of 
high-velocity stars, while the width of the Gaussian part is simi- 
lar. This feature is in agreement with expectations: High-velocity 
stars are usually the lowest mass members of temporary few- 
body systems which are expelled after a close encounter (see 
§3.2\ . The close passage and large mass of the perturber resu lts 
in a significant removal of disc material dPfalzner et al.l l2006). 

In analogy to Fig. [3] positions and velocities of high-velocity 
stars from the numerical simulations are displayed in Fig. [6] 
Here the great advantage of numerical simulations becomes ap- 
parent: several runs of the same model can improve statistics 
far enough to produce prominent features where only weak sig- 
natures in observational data are found. Our dynamical model 
of the ONC reproduces the observed features: As expected, 
most stars are concentrated in the inner tenth parsecs around 
the cluster centre, while several stars are located in the outer 
cluster parts, moving in radial directions from the cluster cen- 
tre. However, we find the same two unexpected features as in the 
observations, namely that (i) the region between the cluster cen- 
tre and the outer cluster parts is underpopulated by high-velocity 
stars, and (ii) the velocity vectors of a fraction of stars in the 
outer cluster parts do not point away from the cluster centre. We 
will refer to high-velocity stars that leave the cluster on a track in 
radial direction from the cluster centre, i.e. with the velocity vec- 
tor pointing away from the cluster core with radius R[ n * 0. 1 pc, 
as "radial escapers", while high-velocity stars that do not match 
this condition will be called "orbital escapers". In the following 
we explain the choice of the terminology and discuss the two 
classes of high-velocity stars. 

In contrast to observations we can trace the history of these 
stars in the simulation; this gives us the opportunity to investi- 
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gate the reason for this strange configuration. As we will see, (i) 
is a consequence of (ii), thus we will discuss that first. In Fig. [7] 
two exemplary tracks of stars in the outer cluster parts are dis- 
played, those of a radial (Fig. [7}}) and an orbital escaper (Fig.|7^). 
It is evident that the phase spaces of the two stars are completely 
different. The radial escaper was originally located close to the 
cluster centre, experienced several encounters in multiple pas- 
sages of the cluster centre and was finally expelled in a close 
three body encounter. After the ejection, it is moving on a track 
in a radial direction from the cluster centre, i.e. with the velocity 
vector pointing away from the cluster core. The orbital escaper 
stems from the outer cluster parts, passed on a non-closed orbit at 
a minimum distance of some tenths of a parsec around the cluster 
centre (without significant encounters) and was accelerated suf- 
ficiently by the central mass to leave the cluster on a hyperbolic 
orbit. Most of the time its velocity vector is not pointing away 
from the cluster core. This is most evident at large distances from 
the cluster centre. Only at the two short periods of cluster centre 
passage, i.e. when the radius vector is approximately normal to 
the velocity vector, the determination of the direction of motion 
is eventually not sufficient to discriminate between a radial and 
orbital escaper. The fact that the orbital escaper is leaving the 
cluster although it was initially energetically (but only weakly) 
bound to the cluster is due to a varying cluster potential on a 
time scale shorter than its revolution around the cluster centre. 2 
The main sources of the potential variation are mass segregation 
and the evaporation of (preferentially low-mass) cluster stars. 
Consequently the orbital escaper is accelerated more strongly 
after the second passage of the cluster centre and is less decel- 
erated in the outer parts due to the lower total cluster mass and 
more extended cluster potential. 

The existence of two different classes of tracks explains why 
the region between the cluster centre and the outer cluster parts 
is devoid of high-velocity stars. Radial escapers leave the cluster 
on a very short time scale, while orbital escapers reside for much 
longer in the vicinity of the cluster centre due to their wide non- 
closed orbits. Hence the disjoint spacial groupings of stars after 
1 Myr of cluster evolution are the consequence of disjoint sets of 
initial phase space volumes. 

We find from our simulations that the distinct dynamics of 
the two exemplary high-velocity stars characterize in general 
the dynamics of radial and orbital escapers: they belong to dy- 
namically distinct groups. Moreover, the two groups as well can 
be separated due to their disc properties: a large fraction of ra- 
dial escapers is discless, while most orbital escapers are star-disc 
systems. This morphological distinction is a consequence of the 
dynamical bisection. The morphological bisection translates ob- 
servationally into a photometrical bisection, i.e. the stars would 
be devided into two groups according to the presence of excess 
emission. Such a trend, though only weak, is also present in the 
observational data. 

In Fig. the three isolated stars with excess emission are 
identified as orbital escapers. The situation is more diffcult for 
the two stars with excess emission immediatly below the clus- 
ter centre. As mentioned above, stars close to the cluster cen- 
tre cannot be uniquely identified as radial or orbital escapers 
from the direction of motion alone. Moreover, we have only 
two-dimensional spacial and velocity information, so the true 
distances to the cluster centre and velocities are not known. If 



2 The time scale of the potential variation is related to the crossing 
time of the cluster, ? cross a 0.5 Myr, while the time scale of revolution 
is approximately the cluster age, ?onc ~ 1 Myr, which is about twice as 
large. 
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Fig. 7. Tracks of two high-velocity stars from simulations. The 
cross marks the cluster centre. The solid line marks the trajec- 
tory up to 1 Myr, the dotted line represents later times, a) Star on 
a wide orbit around the cluster centre, escaping from the clus- 
ter after a close passage of the cluster centre, b) Star escaping 
after multiple passages of the cluster centre and a final close en- 
counter. 



the projected distance and true distance to the cluster centre dif- 
fer only slightly for both stars, then we would classify the more 
distant (~0.25pc) as a probable orbital escaper. The reason is 
that due to its proximity to the cluster centre and relatively low 
velocity this star will be accelerated and deflected by the central 
cluster mass and pass on a curved trajectory. When passing the 
cluster outskirts, its direction of motion would not point away 
from the cluster core and thus it would be identified as an or- 
bital escaper according to our classification scheme. Of course, 
we cannot exclude the possibility that this star was ejected in a 
three-body encounter, although it could not be classified as a ra- 
dial escaper due to its predicted trajectory. For the other star a 
classification as a radial escaper seems more appropriate. If, on 
the contrary, the true distance is much larger than the projected 
distance for both stars, they would be classified as orbital esca- 
pers. However, we can only speculate about the dynamical origin 
of the two stars. The special case of the two close-by stars with 
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excess emission in the lower left will be discussed below. The 
only star without excess emission (JW 45) seems to have been 
expelled very close (<0.13pc) to the cluster centre and is thus 
classified as a radial escaper. If we thus interpret, conversely to 
the previous arguments, the lack of excess emission as an indica- 
tor of the absence of a disc, then this star provides evidence for 
encounter-triggered disc destruction. The signature of disc ma- 
terial of the other stars is - as far as a classification is possible - 
in accordance with our numerical results and thus supports this 
view. 

The two close-by stars at approximately (- 1 .5 pc, - 1 .0 pc) 
in Fig. [3] JW 3 and JW 4, seem to form a binary. Though 
their separation of about 10 4 AU is large, the remarkably sim- 
ilar proper motion s and radial veloci ty (v^ 3 =29.1 kms -1 , 
i/J W4 = 31.6kms _1 : Stassun et al. 1999|), age, and infrared ex- 
cess strongly support the assumption of a physical pair - at least 
in the past. If this pair was expelled as a binary from a four- 
body encounter, than this must have occurred less than 0. 1 Myr 
ago (accounting for the actual distance from the cluster centre, 
the velocity and the deceleration by the interior cluster mass). 
The difference in proper motion of the two stars corresponds to 
a distance of ~ (4 ± 2) • 10 4 AU, in good agreement with the ob- 
served projected separation. Due to the direction of motion, the 
two stars are classified as radial escapers. The expulsion of bina- 
ries from close four-body encounters in our simulations, though 
a rare event (9 events from our 20 runs), usually does not lead 
to a significant disc-mass loss of the individual stars. This again 
is in good agreement with the observational data. Thus the ex- 
cess emission of the binary radial escaper does not contradict 
the correlation of our dynamical and the photometric classifica- 
tion, which in the case of single stars shows that in most cases 
radial escapers are discless and orbital escapers are found pref- 
erentially among star-disc systems. 

The general difference of the orbits of discless stars and star- 
disc systems is represented in Fig. [8] The minimum distance 
of star-disc systems to the cluster centre d m i n is clearly a linear 
function of the initial specific angular momentum /rj (relative to 
the cluster centre). From a linear least squares fit we find a slope 
a with a small standard error Aa, a — (0.902 ± 0.044) s kirT 1 
(Aa/a = 0.048). This means that Iq is conserved - a conse- 
quence of the wide orbit around the cluster centre without strong, 
abrupt perturbations from single stars. On the contrary, disc- 
less stars show a much wider, not clearly correlated distribu- 
tion of dmin due to angular momentum exchange in close en- 
counters in the cluster centre. Here the slope (3 of the linear 
best-fit has a large standard error A/3, /3 = (0.571 + 0.127) s km -1 
(A/3//3 = 0.222). Those star-disc systems with lower Iq and d m \ n , 
populating the discless regime, are components of binaries. 

The distinct dynamics are even more evident from Fig. [9] 
here we show the ratio of the actual (at 1 Myr) and initial specific 
angular momentum l n0 vi/h) as a function of the initial specific 
angular momentum Iq. Star-disc systems are concentrated nearly 
symmetrically around / now //o = 1, while for most discless stars 
Inow/lo > 1, and even up to several tens. The increase of angular 
momentum of the discless stars in relation to the cluster centre 
can be explained as follows: Single stars gain a l arge amount of 
angul ar momentum in a close triple encounter d Valtonen et alj 
2005), and leave on straight radial tracks after breakup. Since 
these encounters occur preferentially close to the cluster centre, 
on average the angular momentum relative to the cluster cen- 
tre is also highly increased. The raise of l now /lo with lower Iq 
is given by the fact that close encounters are more probable for 
stars with lower Iq, which in turn lead to a higher gain in angular 
momentum. 
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Fig. 8. Minimum distance of the high-velocity stars to the clus- 
ter centre, d m i n , as a function of the initial specific angular mo- 
mentum, Iq (relative to the cluster centre). Star-disc systems are 
marked by open, discless stars by filled circles. Linear best-fits 
of both populations, star-disc systems (dotted line) and discless 
stars (solid line), have been included (quantitative results are dis- 
cussed in the text). 
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Fig. 9. Relative actual angular momentum of the high-velocity 
stars, l now /lo, as a function of the initial specific angular mo- 
mentum lo (relative to the cluster centre). Star-disc systems are 
marked by open, discless stars by filled circles. 



The additional components in our numerical model of the 
ONC - a background potential and primordial binaries - have 
different effects on the sample of high-velocity stars. By com- 
paring with results of our simulations without either one or both 
additional components, we find that (i) the inclusion of primor- 
dial binaries has the effect of increasing the number of discless 
stars (due to a larger number of encounters) and to increase the 
maximum velocity of the high-velocity stars (due to a higher 
probability of closer encounters), (ii) a background potential re- 
duces the number of discless stars and reduces the maximum 
velocity of the high-velocity stars, because the higher velocity 
dispersion of the stars reduces the probability of close encoun- 
ters. The combination of both as in the present model does not 
cancel out the positive effect of the binaries and thus results in 
a higher number of discless stars and a higher maximum veloc- 
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ity of the high-velocity stars than would result from a single star 
model. 



5. Summary and discussion 

Combining observational data and numerical simulations, we 
have shown that even after 1 Myr of dynamical evolution a clus- 
ter of the size and density of the ONC is still dynamically active. 
Young stars, preferentially of low-mass, are expelled in close 
/V-body encounters, losing a large fraction of their circumstellar 
matter. 

The encounters are in most cases interactions with mas- 
sive stars in the cluster centre. This finding highlights the sig- 
nificant effect that encounters in (massive) young stellar clus- 
ters can have on the evolution of protoplanetary discs. This is 
even more evident if one addresses not only the effect on the 
mass, but the even stronger effect on the ang ular momentum 
of the disc as shown in p revious investigations (Pfalzner 2006; 
iPfalzner & Olczakll2007al) . 

Using the observational data from I Jones & Walker! d 1988b 
we show (Fig. QJ) that in the ONC there is a small population 
of stars with proper motions // x y larger than three times the one- 
dimensional velocity dispersion of the cluster, fi x y > 7 . 5 km s ~ 1 . 

From numerical simulations we find that the so-called "high- 
velocity stars" form two dynamically disjoint groups. One group 
is composed of stars expelled in a close encounter, moving on 
radial tracks directly outward from the cluster centre, termed 
"radial escapers". The other contains unperturbed stars running 
on wide, non-closed orbits around the cluster centre (on average 
with lower velocities), termed "orbital escapers". The dynam- 
ical state of the stars has its origin in their initial phase space 
location: radial escapers are initially located close to the cluster 
centre (r <0.3pc) and characterized by a low angular momen- 
tum (relative to the cluster centre). In contrast, orbital escapers 
were formed in the outer cluster with large angular momentum 
relative to the cluster centre. 

The different dynamics of the high- velocity stars implies a 
signature in their circumstellar characteristics: stars being sub- 
ject to close encounters are ex pected to lose their d isc material 
faster and to a higher degree dPfalzneretaIi r2006). Using our 
prescription of disc-mass loss in stellar encounters, indeed we 
find a clear correlation from numerical simulations: radial esca- 
pers, initially located close to the cluster centre and later expelled 
in a close encounter, lose more than 90 % of their disc material 
in 1 Myr of dynamical evolution, while orbital escapers, initially 
distant stars, moving on wide orbits around the cluster centre, do 
not experience strong perturbations and keep most of their disc 
material. 

We compare our num e rical results wi t h inf rared ob servations 
fromlHillenbrandl (1 19971) . iMuench et all d2002l) . and lLada et al l 
(12004 7 tracing the (inner) circumstellar material. Though strict 
conclusions are not possible due to the small observational sam- 
ple and some stars that cannot be classified according to our 
scheme, we observe the same trend: stars that have been clas- 
sified as orbital escapers do show near-infrared excess emission, 
indicative of circumstellar matter, while those classified as ra- 
dial escapers show pure photospheric colors, lacking evidence 
of (inner) circumstellar discs. 

This view is strongly supported by the work of iTanl d2004l) 
on the dynamics of 8 l C Ori. The most massive star in the ONC, 
9 l C Ori, has a proper motion several times greater than the dis- 
persion of bright ONC stars and much larger than the velocity 
expected if it were in equipartition with the other cluster stars 



(Ivan Altena et al.lll988l) . lTanl d2004l) showed that the direction of 
8 l C Ori's motion is consistent with being exactly opposite to that 
of a B-type star embedded in the background molecular cloud, 
the so-called BN object, and concludes that it is most probably a 
runaway star originating from the Trapezium ejected about 4000 
years ago after a close encounter with l C Ori. This supports the 
idea that stellar encounters may be likely events in such dense 
regions as the ONC, or at least in their cores. The evidence that 
the most massive stellar object was involved in a close encounter 
with another massive star is in best accordance with the w ork of 
iMoeckel & Ballvl (l2007al) and lPfalzner & Olczakl(l2007bl) . 

Considering the dynamical age of the ONC of several cross- 
ing times, encounters must have had much stronger impacts on 
a stellar disc at an early age of the cluster, when densities were 
much higher but massive stars already had formed. Indeed, our 
simulations confirm this expectation, giving rise to an era of 
strong and frequent interactions among star-disc systems at the 
onset of massive star formation. 

We expect that in even denser clusters such as the Arches 
cluster, high-velocity stars should be even more frequent than in 
the ONC. In such systems the current spatial and velocity distri- 
bution of the high-velocity stars should give strong indications 
of how the cluster developed in former times. 
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Appendix A: Observability of tidal tails due to 
star-disc encounters in the ONC 

The relevant physical quantity that determines the prominence 
of tidal tails due to an encounter i s the change of angular mo - 
mentum in the disc. According to iPfalzner & Olczakl d2007af) . 
we assume that a fractional angular momentum loss (AML) of 
at least 10 % is required to for m observationally detec table tidal 
tails (see also Fig. 9 and 10 of Pfa lzner & Olczakll2007al) . 

To have an estimate of the rate of encounters in the ONC in 
which the AML is at least 10 %, we consider an encounter of a 
star with mass m = 0.5 M , which corresponds to the mean stel- 
lar mass in the ONC (see A star of that mass is thought to be 
surrounded by a plotoplanetary disc of typical size r& = 100 AU. 
To be on the safe side, we want to estimate the upper limit of 
the encounter rate and further assume the encounter partner to 
be the hig hest mass star of th e ONC, M = 50M o . Referring to 
Table 1 of IPfalzner & Olczakl d2007al) . the specified minimum 
AML requires an encounter at a relative periastron r p /rd~ 10, 
or a periastron r p a 1000 AU. Moreover, we assume that the en- 
counter occurred in the Trapezium Cluster (TC), the dense cen- 
tral part of the ONC , where it is most probable (see Fig. 4 of 
Pfalz ner et alj|2006b . The number of stars, the density and the 
velocity dispersion of the TC are Njc ~ 750, «tc ~ 10 3 pc~ 3 , and 
cr j w = 2.5kms- 1 (see ©. 



The time scale for encounters of the assumed type is given by 
Eq. dC.ll i; substituting n = rtjc, c = °"id> r * ~ r p> anc ^ m * ~ M, we 
arrive at f enc « 2 • 10 5 yr. With a dissipation time scale of the tidal 
tails of fdi ss < lOOOyr, the probability of a detection is roughly 
^obs 5 fdissAenc ~ 5 • 10 3 . With the number of stars in the TC, 
A^tCj we expect at most N b s ^ PobsNjc ~ 4 stars to be accom- 
panied by tidal tails that could be observed at the current time. 



Appendix B: Estimate of the mean uncertainty of 
stellar ages 

Since ind ividual errors of th e derived stellar ages are not pro- 
vided bv lHillenbrandldl997l) . we estimate the mean error from 
the quoted observational and theoretical uncertainties and the 
constructed HR diagram. The uncertainties of the derived lumi- 
nosities, log(L»/L Q ) < 0.2, translate into an age uncertainty of 
~0. 2-0.4 dex for stars with masses 1.0-0.1 M . Uncertainties of 
the derived effective temperatures, log T e g < 0.02, translate into 
an age uncertainty of typically ~0.3 dex, but can be as large as 
~1 dex for stars with mass <0. 1 M©. Further uncertainties of the 
derived ages are introduced due to differences between different 
pre-main sequence evolutionary tracks, which can be as large as 
0.6 dex. Assuming a typical uncertainty of 0.3 dex due to uncer- 
tainties from luminosity, effective temperature and evolutionary 
tracks we estimate a mean uncertainty of stellar ages of ~0.5 dex. 



Appendix C: Minimum velocity for unperturbed 
escape of stars in the ONC 

High-velocity stars that have been generated due to a close triple 
encounter in the cluster centre are expected to leave the clus- 
ter without significant perturbation. This is due to the fact that 
(a) the fractional change of the velocity u of a high-velocity star 
is less t han 10% unless the impac t parameter is not lower than 
100 AU dBinnev & Tremaind 19871 Eq. 4-8), and (b) only a small 
fraction of stars experiences more than on e encounter closer than 
lOOAUin one cros sing time of the ONC dScallv & Clarkd2001l: 
Olczak et al. 2006). Alternatively, one can evaluate the collision 
time scale of the escap er and a binary with se mi-major axis 
a w 100 AU in the ONC dBinnev & Tremaind[l987l Eq. 8-123), 



*coll 



16 ^fnncrrl ( 1 + 



2o- 2 r 



lOMyr, 



(CI) 



where the gravitational focusing of a 1OM binary in the dense 
Trapezium Cluster has been considered in order to obtain a ro- 
bust lower limit on f co n. We have used r t -a,n- n T c ~ 10 3 pc~ 3 
and cr = <t^ w = V30 -1 * = 4.3 km s _1 (see ©. Since the collision 
time scale is much longer than the time to reach the cluster out- 
skirts, foot] > lOMyr » 0.2 Myr > f eS c, it is valid to assume that 
stars with velocities v > 3o-^* are effectively unperturbed before 
they escape from the cluster. 



Appendix D: Maximum velocity of a star ejected 
from a bound triple system 

We analyze a triple system after the ejection event forming a con- 
figura tion of a binary and an escaping body (cf. IValtonen et al.l 
120051) . The total energy of the system is then 



E = -mr - G-^ + -Mr 2 - G—^—, 
2 s r s 2 r 



(D.l) 
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where r s is the separation of the third body relative to the 
barycentre of the binary, r the separation of the binary compo- 
nents, m s , m- d and m\, the masses of the escaper and the binary 
components, = m. d + m\, the binary mass, and M = m a mb/#iB 
and m = m^m % l '(ma + nh) the reduced masses. 

When observing, the detection of the escaper will usually oc- 
cur when the distance to the binary components greatly exceeds 
the size of the binary system, r s » r. Thus we can neglect the 
second term in Eq. ( ID. lb . Moreover, in most cases the ejected 
body will be the lowest mass component of the triple system, 
m s <sc otb (see 33.21 1. hence Eq. ( ID. 1ft reduces to 

1 i 1 n trinTnh 1 m a mi, 

E = -m s r 2 + -Mr 2 - G—^—^- = -m s f 2 - G^-A (D.2) 

2 2 r 2 2a 

where the last term denotes the total binary energy in the general 
case of an elliptical orbit with semi-major axis a. 

If we use the same parameters for the triple system as 
in 33.21 and assume a nearly equal-mass binary, m s = 0.25 M , 
2a — 100 AU, m a mb~4M , we find from Eq. dP.2| i the maxi- 
mum velocity of the escaper from a triple with negative total 
energy, v s — r s < 18kms _1 . In fact, our simulations show that 
the binary that generates a high-velocity star never is less mas- 
sive then 10 M and usually exceeds 20 M Q , and the mass ratio is 
never below 1/6 and usually about 1/3. Using these parameters, 
we find v s < 30-80 km s" 1 . 
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